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Abstract 

We propose Raman-induced collinear difference-frequency generation (DFG) as a method to 
manipulate dynamical magnetization. When a fundamental beam propagates along a threefold 
rotational axis, this coherent second-order optical process is permitted by angular momentum 
conservation through the rotational analogue of the Umklapp process. As a demonstration, we 
experimentally obtained polarization properties of collinear magnetic DFG along a [111] axis of a 
single crystal of antiferromagnetic NiO with micro multidomain structure, which excellently agreed 
with the theoretical prediction. 

PACS numbers: 78.20.Ls, 42.65.-k, 76.50.+g 
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Recent developments in pulsed laser techniques have enabled ultrafast optical manipu- 
lation of magnetization in solids [l-lOl. In particular, coherent control through Raman- 
induced nonHnear optica, processes ffl has attracted interests, because low- euergy mag- 
netic excitations are accessible without excess heating |5J. Therefore, compared with the 
conventional processes involving optical absorption followed by thermal relaxation 6|-[9), 
these coherent processes have potential applications including quantum information pro- 
cessing and ultrafast switching techniques for spintronics. 

In such coherent optical processes, conservation of angular momentum determines the 
polarization of the induced excitation, whereas conservation of energy and momentum de- 
termines its frequency and wave vector, respectively 
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12] . The balance of angular mo- 



mentum includes contributions from the electromagnetic field (J EM ), excitations in solids 
(J cx ), and the crystalline lattice (J c ): AJ EM + AJ CX + AJ C = [12]. When the z-axis is a 
continuous rotational axis, J EM + J| x must be rigorously conserved because J c z is a conserved 
quantity under any rotation around the z-axis. By contrast, along a threefold rotational axis 
of the crystal, J EM + J| x is conserved within a change of 3H] this is the rotational analogue 
of the Umklapp process 12j, ll3[. This discrete rotational symmetry in solids, and the resul- 
tant polarization selection rules in nonlinear optics, has been taken into account for decades 
when the second-harmonic generation [12I . ["K^ and the parametric down-conversion [14] are 
discussed. However, little attention has been paid to implications of the discrete rotational 
symmetry for the coherent control of elementary excitations in solids. 

In this Letter, we discuss selection rules for optically induced magnetization from a view- 
point of conservation of angular momentum. The inverse Faraday effect 15] and the mag- 
netic difference-frequency generation (DFG) are used as specific examples. Both are de- 
scribed by third-rank axial c-tensors ^( 2 ) MEE that change sign under time reversal [if], \v\ . 
The superscripts label magnetic (M) or electronic (E) interactions. Subsequently, we pro- 
pose collinear magnetic DFG along a threefold rotational scheme of optical control 
of magnetization, using the discrete rotational symmetry. As a demonstration, polarization 
selection rules are calculated and experimentally obtained for a collinear DFG in NiO that 
involves two frequency components of an ultrashort laser pulse, showing excellent agreement. 

To discuss the conservation of angular momentum in light-matter interactions for photons 
propagating along the z-axis, it is convenient to distinguish vector fields by their helicity. 
Any electric field E with a propagating vector k along the z-axis can be decomposed into 
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circular bases as: E(u;, z) = 2Re [(£' + e + + £Le_ + E z e z ) e -nwf-ft*)J ) with unit vectors, e + = 
(e x +ie y ) j V2, e_ = (e x — ie y )/\/2, and e z . E + and E_ are left- and right- circular light fields, 
whereas E z is longitudinal. As illustrated in Fig. [U^a), creation or annihilation of a photon 
induces A J™ = ±h in correspondence with its helicity, while Jf M = for the longitudinal 
field along the z-axis. The same rule is available for a magnetic field represented in circular 
bases. Therefore, if and only if the angular momentum of an elementary excitation is ±h, 
its transition to and from one photon that propagates along the z-axis is allowed. 

A light beam at the frequency u can induce dc magnetization M(0) through the inverse 
Faraday effect: 

Mf >(0) = xgi MEE (0; -u,u)E*(u)E k (uj), (1) 
where E(w) is the electric field of the light wave [15] . In a nonabsorbing isotropic medium 

(2) 

invariant under time- reversal, x\j k = X £ ijk- Here, x ls a scalar and e^k is the Levi-Civita 
tensor. Then, Eq. (P) is reduced to M^ 2 ) = %E* x E = xd-^+l 2 — l^-l 2 )^. Thus, in an 
isotropic medium, an elliptically polarized light beam induces M^ 2 ^ along its propagating 
direction. 

Similarly, ac magnetization M( 2 )(£7) at the difference frequency f2 = uj x — uj 2 can be 



parametrically generated by two lig 
Raman-type nonlinear process jl8|, 



it waves with frequencies U\ and u 2 through a stimulated 
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M? ] (Sl) = x%l MEE (n; -u u u 2 )E*ME k {u 2 ). (2) 

When M^ 2 ) is not only Raman (two-photon) active but also infrared (one-photon) active 
(via a magnetic dipole process), magnetic field H^ 2 ^(f2) and electric field E^ 2 ^(^) are radiated 



from M^iyt) [18]: 



Q 2 

V x (Vx) - — e r 



c 
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H (2 \n) = — e r M (2) (fi) (3) 



c 2 



-ine r e E {2) (n) = V x H^(fi), (4) 

where e r is the relative permittivity at Q, e Q is the vacuum permittivity, and c is the speed 
of light; this is the magnetic DFG. 

In continuous rotational symmetry, as in a liquid or gas, the requirement of rigorous 
conservation of Jf + J| x forbids collinear DFG — i.e., DFG with incident and out-going 
photons propagating along the same (z-) direction. In particular, optically induced magnetic 
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excitation in the medium has J| x of either or ±2h, and is thus nonradiative [Figs. QJb) 
and (c)]. Note that generation of longitudinal field (J° x = 0), rather than DFG, is possible 
by scattering the incident photon to a mode with the same helicity. If the scattered mode 
is the identical one with the same frequency, the result is a dc field along the z-axis, which 
contributes to the inverse Faraday effect. 

However, DFG is not necessarily forbidden if a fundamental light beam propagates along a 
crystal's threefold symmetry axis and involves two photons of opposite helicity — a linearly 
polarized pulse satisfies this condition, while purely circularly polarized pulse does not. 



This is because this discrete symmetry allows the change of J c z by 3h |12l. |14|]. For example, 
annihilation of the right-circular photon and creation of the left-circular photon induces a 
magnetic excitation with J° x = —2h. This is equivalent to +H through decrease of J c z by 
3h, and emits a left-circular photon [Fig. QJd)]. The result is that Jf M increases by 3h and 
the total angular momentum is conserved. For an incident pulse that is linearly polarized, 
a mirror process occurs simultaneously, where AJf M = —3h and AJ Z = +3h. 

Nonvanishing transverse components of ^-( 2 ) MEE invariant under rotation by 2ir/3 around 
z-axis describes the allowed collinear DFG process along a threefold axis: 

Xxxx Xxyy Xyxy Xyyx — OC, (5) 

Xyyy Xyxx Xxyx Xxxy = &i (6) 

and other components vanish. Here, a and /3 may also be zero depending on additional 
spatial or temporal symmetry operations allowed. In the circular basis, they are 

X+++ = X^(e + )*(e + )*(e + )* fe = V2{a + (7) 
X— = Xy *(e_)?(e_);(a_); = V2(a - tf), (8) 

which correspond to the processes that involve change of J™ by 3H and —3ft, respectively. 
When the amplitudes of these two processes are the same, as in the case of NiO discussed 
below, the resultant radiation is linearly polarized for an incident pulse with any fixed 
polarization. 

In the remaining paragraphs, the above description of the collinear DFG will be applied to 
the interpretation of our experimental observation of THz radiation from antiferromagnetic 
magnons in NiO illuminated by laser pulses. We chose the direction of the propagation of 
light to be parallel to the [111] axis of an as-grown crystal of NiO that consists of randomly 



distributed domains due to magnetic ordering. This axis is effectively threefold if the signal 
induced in each domain is coherently superposed, as we will discuss below. NiO has a 



magnon resonance around 1 THz that is both Raman [20( and infrared [2l|, |22J active. 
Recently, researchers reported radiation from the magnon [3] and oscillation of the Faraday 
rotation |4| after illumination by pulsed lasers. Our setup differs from these previous studies 
in the selection of the crystal orientation and the domain structure because it is crucial 
for a clear demonstration of this DFG process. Note also that the absence of spontaneous 
magnetization in NiO rules out photomagnetic effect [l Of ] . 

NiO crystallizes in a rocksalt structure above its Neel temperature of Tn = 523 K. Below 
T N , the Ni 2+ ions with opposite sign of spins (S = 1) occupy equivalent crystallographic 
positions [Fig. 0a)], forming two sublattices. Subsequent lattice distortion induces four 
variants called T-domains, and each of them have three variants called S-domains that 
differ by their spin orientation (total of 12 variants) 23] . For one of the T-domains, T x , the 
crystal shows slight rhombohedral distortion along its [111] axis; for T2 to T4, the distortion 
is along the other three diagonals. The domains are accommodated by twin structures 
randomly distributed in an as-grown NiO crystal [Fig. EJ^b)] 24j, and their sizes were several 
micrometers according to polarization microscopy. The spin orientation in one of the three 
S'-domains belonging to Ti-domains (denoted by Ti(Si)) is [112] and for Ti(S < 2) and Ti(S^), 
their equivalents in the (111) plane [Fig. EJa)]. 

To examine whether all domains coherently contribute to the magnetic response, we 
measured the magnetic susceptibility tensor of an as-grown NiO crystal near its magnon 
resonance (~ 1 THz) by time-domain spectroscopy at room temperature [Fig. EJ^c)]. Ellip- 
sometric measurements showed that the magnetic response is isotropic [25] . This indicates 
that although the susceptibility of each 5-domain is anisotropic Q , coherent superpo- 
sition of the 12 variants results in effectively isotropic susceptibility, described by a scalar 

The magnetic susceptibility of this magnon resonance is calculated as: X^) — 
1 7 2 H A {2AM+Hj!^-u) 2 +2iTu> ' wnere 7 is the gyromagnetic ratio of an electron, V is the damp- 
ing coefficient, M is the sublattice magnetization, and AM and are the exchange and 
anisotropy fields, respectively 0, The factor of 2/3 is the result of averaging over the 
domains. Figure EJ^d) is the plot of x^i 00 ) using M, A and if a from the literature Q] and 
taking thermal suppression of M into account [2l|, |22J. The quantitative agreement with 
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the experimental result assured us that the S'-domains are sufficiently large to accommodate 
magnons. 

The DFG in this sample, therefore, is also described by coherent superposition of the 
signal generated in each of the 12 variants: 

n 

Here, -^-( 2 ) MEE > n i s the nonlinear susceptibility tensor of an n-th domain, which can be char- 
acterized by its orientation with respect to the crystal's Cartesian frame, denoted by su- 
perscript (C), as depicted in Fig. [3]^a). In a single T 1 (S'i)-domain, the allowed symmetry 
operations are 1, 1, 2 y , and 2 y , and magnetic point group of a Ni site is 2/m (the underbars 
denote time- reversal operation) This gives 14 nonvanishing independent components 

of x {2)MEE,n , including Xzzz, Xxxx, Xyyx, Xzzx, Xxxz, Xyyz, and other components obtained by 
permutation of the Cartesian indices. Tensors for other variants are obtained, as is discussed 
for magnetic second-harmonic generation |24|: the tensor of a Ti^) -domain is calculated 
by transforming ^( 2 ) MEE A under rotation by 2%/3 around the z( c )-axis. The nonvanishing 
transverse components of ^( 2 ) MEE in the crystal coordinate frame is derived by averaging 
^,(2)MEE,n over ^ variants as: 

Xxxx = Xxyy = Xyxy = Xyyx = Oi. (10) 

These are precisely the transverse components of a tensor that is invariant under rotation 
by 27r/3: a special case of Eqs. (jSJ) and ([6]) in which /3 = 0. Therefore, this multidomain 
crystal behaves in the same way as a single-domain crystal with threefold symmetry as far 
as the DFG is concerned, and thus collinear DFG along the z^-axis is allowed. 

We performed time-domain spectroscopy measurements of THz radiation from NiO. The 
crystal was irradiated by Ti-sapphire-based ultrashort laser pulses (~ 100 fs) along its [111] 
crystallographic shown in Fig. [3]Ja). With linearly polarized excitation, the DFG was 

observed, but not with circular polarization. Figures Mjo) and (c) show a typical temporal 
waveform of the generated THz field E( 2 )(t), and its Fourier transformed spectrum |E^ 2 ^(f2)|, 
respectively. Radiation with a resonant frequency of Qq ~ 1 THz was observed. Note that 
;he absorption coefficient of NiO is small (60 cm -1 ) for the laser wavelength (~ 800 nm) 



29(, and thus, the theory for nonabsorbing media is applicable. 

Polarization properties of the THz radiation were observed by THz ellipsometry. Incident 
light E p was fixed to be x^-polarized, and the sample was rotated by an angle 9 around the 



z( L )-axis: the superscript (L) denotes the laboratory coordinate, as shown in Fig. [3]^a). The 
radiation was linearly polarized for any 9, and the x^- and y^- polarized components of 



E^ 2 ) (t) were acquired using a wire-grid polarizer installed behind the sample 25j . Transverse 
components of M( 2 )(f2 ) at the resonant frequency fio were determined as functions of 9, 
using Eqs. and (jl]), and are shown in Figs. [3]^d, e). 

In the laboratory coordinate frame, M x and M y can be calculated by Eq. (j2j), in 
which the susceptibility tensor is derived by rotating the tensor in Eq. ( TTU1) by 9 around the 
z^-axis: 

= a|E p | 2 cos3fl, = a|E p | 2 sin3#. (11) 

Here, E p || ei L ^ is used. Angular dependence of Eq. ffTTj) is depicted in Figs. [3](d, e), showing 
excellent agreement with the experimental data. In addition, the measured amplitude |M^| 
was proportional to |E p | 2 , indicating that the radiation mechanism is indeed a second-order 
nonlinear optical process. It is important that the THz ellipsometry along a linearly isotropic 
axis allows us to obtain directly the orientation of the induced magnetization in contrast 



with nonlinear optical techniques 



r optK 

0,9, 



the 



30]; for example by the Faraday rotation 
Kerr rotation [5|, IfSj, or second-harmonic generation 9]. 

In summary, we showed that magnetization can be induced through a collinear magnetic 
DFG process, which is allowed with linearly polarized excitation along a threefold axis. 
The rotational analogue of the Umklapp process plays a key role because it allows the 
change of Jf M + J| x by ±3h. We experimentally demonstrated the DFG resonant on the 
antiferromagnetic magnons in a multidomain NiO crystal. We observed clear polarization 
dependence on the incident and radiated photons under linearly polarized excitation, in 
excellent agreement with the theoretical description. It has highlighted that the conservation 
of angular momentum can provide a general guiding principle to design experiments by 
taking material's symmetries into account. Note that this concept is readily applicable to 
other elementary excitations induced through the nonlinear optical processes, for example, 
electronic DFG or optical rectification because both electric and magnetic fields behave 
in the same way under a rotation that maintains the parity (proper rotation). Extensive 
applications of this approach to other systems, such as artificial structures or atomic gasses 
in optical lattices, are expected. 
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FIG. 1: (Color online) (a) Changes of J^ M by creation or annihilation of photons propagating 
along the z-axis. The lower panels show schematics of the changes of the angular momentum in 
collinear scattering processes of photons to a mode with (b) the same helicity and (c) the opposite, 
and (d) a collinear DFG process allowed along a threefold axis involving two photons with opposite 
helicity. 
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FIG. 2: (Color online) (a) Schematics of spin alignment in three S'-domains belonging to the T\- 
domain, and (b) random distribution of four orientations of T-domains within the coherent length 
Ac of the light beam, (c) Experimental and (d) calculated magnetic susceptibility of NiO; T = 70 
GHz and e r = 12.25 are assumed. 
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FIG. 3: (Color online) (a) Schematics of the experimental setup. Inset shows the relation between 
crystal and laboratory coordinates, (b) Temporal waveform of the radiation from NiO and (c) its 
Fourier transformed spectrum. Polar plots of the transverse components (d) M x and (e) M y of 
the induced magnetization as functions of the crystal rotation angle 9. The dots are experimental 
data of their amplitude, and the solid curves are calculated by Eq. (fTTjh 
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